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ABSTRACT 

Black hole feedback - the strong interaction between the energy output of su- 
permassive black holes and their surrounding environments - is routinely invoked 
to explain the absence of overly luminous galaxies, the black hole vs. bulge cor- 
relations and the similarity of black hole accretion and star formation histories. 
Yet direct probes of this process in action are scarce and limited to small samples 
of active nuclei. In this paper we present Gemini Integral Field Unit observations 
of the distribution of ionized gas around luminous, obscured, radio-quiet quasars 
at z ~ 0.5. We detect extended ionized gas nebulae via [O lll]A5007A emission m 
every case, with a mean diameter of 28 kpc. These nebulae are nearly perfectly 
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round, with H/3 surface brightness dechning oc R~^'^^^'^. The regular morpholo- 
gies of nebulae around radio-quiet quasars are in striking contrast with lumpy or 
elongated [O ill] nebulae seen around radio galaxies at low and high redshifts. 
We present the uniformly measured size-luminosity relationship of [O ill] nebulae 
around Seyfert 2 galaxies and type 2 quasars spanning six orders of magnitude in 
luminosity and confirm the flat slope of the correlation (i?[o m] oc Lj'q^^j^^ ^^). We 
propose a model of clumpy nebulae in which clouds that produce line emission 
transition from being ionization-bounded at small distances from the quasar to 
being matter-bounded in the outer parts of the nebula. The model - which has 
a declining pressure profile - qualitatively explains line ratio profiles and surface 
brightness profiles seen in our sample. It is striking that we see such smooth and 
round large-scale gas nebulosities in this sample, which are inconsistent with il- 
luminated merger debris and which we suggest may be the signature of accretion 
energy from the nucleus reaching gas at large scales. 

Subject headings: quasars: emission lines 



Introduction 



Feedb ack from accreting black holes has beco n ae a key element in modeling galaxy 
evolution (ITabor &: Binneyl Il993l : ISilk &: Reed Il998l : ISpringel et al.l l2005l ) . For one thing, 
regulatory mechanisms are needed to keep massive galaxies from forming too many stars 



and becoming overly massive or blue at late times (IThoul fc Weinberg] 1 19951 : ICroton et al. 
20061 ). Black holes are natural candidates because even a small fraction of the binding 



energy of the infalling material is in principle sufficient to liberate the galaxy-scale gas 
from the galaxy potential. The discovery that the masses of supermassive black holes 
in inactive galaxies strongly correlate with the velocity dispersions and masses of their 



hosts' stellar bulges (IMagorrian et al 



Tremaine et al. I I2OO2I : 



1998; Gebhardt et a 



Marconi fc Huntll2003l : iHaring fc Rix 



2000 



2004 



Ferrarese &: Merritt 



Giiltekin et al 



2009 



2000 



McConnell et al. 



20 111 ) also suggests that the energy output of the black hole in its m ost active - quasar - 
phase must be somehow coupled to the gas from which the stars form (IHopkins et al.ll2006l ). 
The problem remains to find direct observational evidence of black hole-galaxy self- regulation 
and to obtain measurements of feedback energetics. 

In some situations, accretion energy clearly has had an impact on the large-scale en- 
vironment of the accreting black hole. For instance, radio jet activity in massive elliptical 
gala xies and brightest cluster galaxies deposits energ y into the hot gas envelope of the clus- 
ter (IMcNamara &: Nulseru 120071 : iRandall et al.ll201ll ). Likewise, there is clear evidence that 
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powerful radio jets entrain warm 



host galaxies ( 


van Breueel et al.l 


1986; 


Tadhunter 


1991 


1999; 


Nesvadba et al. 


2006[ 


2008; 


Fu & Stockton 


2009|), 



gas a r id carry significant amounts of niaterial out of their 



Whittle 1992; Villar-Martm et al. 



and the observed entropy profiles of 



galaxy clusters can be explained by invoking a heating mechanism from their central radio- 
loud quasars to reduce the gas supply from the hot intra-cluste r medium (IScannapieco fc Oh 



2004 



VoitI l2005l : IVoit fc Donahud l2005l : iDonahue et al.ll2006l ). Nevertheless, as radio-loud 
quasars only represent a small fraction (~ 10%) of the entire quasar population, invoking 
this mechanism as the primary mode of limiting galaxy mass would require all galaxies to 
have undergone a radio-loud phase - a conjecture which lacks direct evidence and contradicts 
a theoretical paradigm in which radio-loudness is determined by the spin of the black hole 
JXchekhovskov et al.ll201o[ ). 



1995 



Nuclear activity is kno wn to drive outflows on small scales ( Crenshaw et al. 20031: 
Ganguly fc Brotherton 20081 ) , perhaps via the radiation pressure of the quasar (IMurray et al. 



Proga et al.ll2000f) . Broa d absorption- line troughs are seen in ~ 10 — 20% of luminous 



quasars ( jReichard et al.l 120031 ) . or alternative ly the outflows are ubiquitous, and this fre- 
quency sirnply r eflects their covering factors ( IWeymann et al.l Il98ll : iGallagher et al.l 12007 
Shen et al.l l2008l ). The velocities in broad absorption lines can be up to ~ 10,000 km/s, 
suggesting that high energies are involved, but in practice velocity is the only physical 
parameter in such outflows which i s directly measur able. Careful photo-ionization model- 
ing of specific absorption features ( lArav et al.l l2008l ) demonstrated that in several objects 
the outflow appears to extend out to se veral kpc from the nucleus and t o carry a large 
amount of kinetic energy and momentum (IMoe et al.ll2009l : iDunn et al.ll2010l ). Such analysis 
requires high-quality ultra-violet spectra of special broad absorption-line quasars showing 
weak metastable transitions and is only possible for a handful of sources. 

Over the past few years, we have undertaken a spectroscopic campaign to use the spatial 
distribution and kinema tics of ionized gas to sear ch for any impact of the black hole on the 
galaxies on large scales ( IGreene et al.ll201ll |2012| ). Our approach is to measure the spatial 
distribution and kinematics of ionized gas in quasars. While such work has b een conducted 
by sev eral groups since it was pioneered by IStockton &: MacKentyl (Il987l ) and iBoroson et al. 
( 119851 ) ■ our approach is different in a couple of respects. First, we have pushed these obser- 
vations to the highest luminosity quasars, because the efficiency of feedback inferred f rom 
simulations of galaxy formation rises steeply with black hole mass ( ICroton et al.ll2006l . al- 
though the most massive black holes are barely active, at least at low redshifts). Second, 
we focused our observations on obscured quasars in which circumnuclear material blocks the 
line of sight to the observer. 



Such objects - the luminous analogs of Seyfert 2 galaxies in the framework of the unifica- 
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tion model (jAntonuccilll993[) - remain ed unidentified in large numbers until the Sloan Digital 



Sky Survey (SDSS; York et al. 200o[ ). Now we have a sample at 2; < 0.8 s elected based on 



emiss ion line ratios and the [O lll]A5007A (h ereafter [O ill]) lin e luminosity ( Zakamska et al. 
2OO3I ) which compris es nearly 1000 o bjects (iReyes et al.ll2008[ ) with bolometric luminosities 
up to 10^'' erg s~^ ( jLiu et al.l 12009). Extensive follow-u p by our group and o thers using 
the Hubble Space 



20061: IVignali et al. 



( IZakamska et al. 



elescope (HST 



2010: Jia et a 



2005h. Gemini 



Zakamska et al.l 120061). Chandra an d XMM ( iPtak et al. 



20121). Svitzer (IZakamska et al.ll2008|). suectropq 



Liu et al.ll2009[ ). V LT (IVillar-Martm et al. 



Alto 3.5m ( Humphrey et al. 201ol ). and the VLA ( Zakamska et al 



2004 



2011 



a 



arimetry 
al)j Calar 



Lai fc Hoi I2OIOI ) 



yields a detailed description of the optically selected obscured quasar population and their 
host galaxies. Among other findings, at the highest luminositie s, it has become c lear that 
these objects constitute at least half of the quasar population (IReyes et al.ll2008[ ) and are 
therefore representative of "typical" black hole activity. These objects allow us to observe 
extended emission line regions without the overwhelming glare of the quasar itself. There 
are even some t heoretical suggestion s that quasars experience outflows preferentially in the 
obscured phase (IHopkins et al.ll2006l ). 



In this paper, we present Gemini-North Multi-Object Spectrograph (GMOS-N) Integral 
Field Unit (IFU) observations of a sample of fourteen quasars and the analysis of the spatial 
extents, morphologies, and physical conditions of their narrow emission line regions. In 
Section [2] we describe sample selection, observations, data reduction and cahbrations. In 
Section El we present maps of the ionized gas emission, in Section H] we discuss physical 
conditions and morphologies of the nebulae and we summarize in Section |3 We use a 
/i=0.71, fim=0.27, Qa=0.73 cosmology throughout this paper. Objects are identified as SDSS 
Jhhmmss.ss-|-ddmmss.s in Table [1] and are shortened to SDSS Jhhmm-|-ddmm elsewhere. 
The rest-frame wavelengths of the emission lines are given in air. 



2. Data and measurements 



2.1. Sample selection 



We study eleven radio-quiet type 2 quasars selected from the catalog by iReyes et al. 



( I2OO8I ) according to the following criteria: 



1. We select the most luminous quasars in the catalog, with [O ill] line l uminosities 
L[o „ r] > 10^^'^ erg s~^, corresponding to intrinsic luminosities Mb < —26.2 ( jReyes et al 



20081 1. 
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2. Among these luminous sources, we select objects with the lowest redshifts {z = 0.4-0.6) 
in order to maximize the spatial information of our observations. 



3. We req uire the radio flux at 1.4 GHz L 



survey ([Becker et al. 



995; White et al 



< 10 mJy as determined by the FIRST 
19971 ). Applying the i^'-correction formula in 



Zakamska et al.l ( 12004J ). this corresponds to uL^, < 1.1 x 10^^ erg s ^ at the 1.4 GHz 



rest-frame frequency for z = 0.5 sources. This is a conservative threshold which puts 
our targets about one order of magnitude below the separation of the radio-loud and 



radio-quiet objects at these [O ill] luminosities and at these redshifts ( IXu et al.lll999 



Zakamska et al.ll2004j ). 



We use the distribution of objects in the [O lll]-Lradio plane for defining radio-loud, 
radio-intermediate and radio-quiet objects because unlike definiti ons based on co ntinuum 
flux, the one based on [O ill] is less sensitive to the optical type. IXu et al.l (119991 ) adopt a 
separation line with Lradio L'^q m] , and follo wing this definition we find that approximately 
10% of type 2 quasars quahfy as radio- loud (IZakamska et al.ll2004j ). although our dynamic 
range is not sufficient to determine whether this fraction varies significantly with luminosity. 
Our adopted d efinition produces similar fractions of radio-loud sources among type 1 and 
type 2 quasars (IReyes et al.ll2008[ ). but yields a lower radio-loud fraction than th e traditional 



defin itions based on a constant optical continuum to radio continuum ratio (e.g. 
2007h . 



Jiang et al. 



In addition to the 11 radio-quiet objects, we observed three radio- loud type 2 quasars 
to provide a comparison sample for ionized gas morpholog function of jet activity. 

Two radio-l o ud ta rgets, SDSS J0807+4946 and SDSS JllOl+4004, were selected from the 
Reyes et al.l (120081 ) catalog according to criteria 1 a nd 2. The FIRST imag e centered on 
the position of SDSS J0807+4946 shows an FR II ( iFanaroff fc RilevI Il974j ) double-lobed 
radio galaxy oriented at 60° East of North. The maximal lobe distance from the host 
galaxy is ~50" (330 kpc). While there is no apparent radio core component, the two lobes 
are symmetric around the SDSS position, and there is no other compelling candidate host 
galaxy for the radio source. Therefore, we identify SDSS J0807+4946 with the radio source 
87G B 0804-1-4955, whi ch has a spectral index of a = —1.0 {F^, oc i/") between 80 cm and 6 
cm (IBecker et al.lll99ll ). 



The radio fiux of SDSS JllOl-l-4004 places it in the radio-intermediate regime (jXu et al. 



19991 : IZakamska et al.ll2004j ). but the FIRST image clearly shows a pair of radio lobes oriented 
at 112° East of North and symmetrically located at ~40" (23 kpc) away from the nucleus, 
which is a point source with a peak fiux of 17.48 mJy/beam (IBecker et al.lll995l ). Using the 
FIRST data, we find the integrated fiux to be 45 mJy and 21 mJy for the northwest and 
southeast radio lobes, respectively. 
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The third radio-loud target is 3C67 at z=0.311, a compact steep-spectrum source with a 
total extent of ~ 3^^ and wit h lobes oriented at 173° East of North in the 5.0 GHz radio band 
( lEracleous fc Halpernll2004f). The optica l spectrum of this object is consistent with type 2 
quasar classification ( ISpinelli et al.ll2006[ ). and its [O ill] line luminosity is slightly lower than 
the majority of the other objects (see Tabled]). The HST/STIS long-sht spectroscopy find s 
the [O III] line emission to extend out to 2.7 kpc away from the nucleus (lO'Dea et al.ll2002l ). 
Not enough information is available for a direct comparison between the sensitivity of these 
observations and ours. We find general agreement between the orientation of the extended 
emission, while our detection of the faint [O ill] emission extends further out, presumably 
because of better sensitivity. 



2.2. Observations and data reduction 

We observe d eleven radio-quiet obscured quasars and three radio galaxies (Table [T]) with 



GMOS-N IFU flAUington-Smith et al.ll2002f l in December 2010 (program ID: GN-2010B-C- 
10, PI: Zakamska). We use the two-slit mode that covers a 5"x7" field of view, translating to 
a physical scale of 30x42 kpc^ at z = 0.5, the typical redshift of our objects. The science field 
of view is sampled by 1000 contiguous 0.2"-diameter hexagonal lenslets, and simultaneous 
sky observations are obtained by 500 lenslets located ~1' away. In view of their respective 
redshifts, 13 objects were observed in the z-band (7060-8500 A) so as to cover the rest 
frame wavelengths 4100 < A < 5200 A and thus H/3 and [O ill]. Of those, two were also 
observed in the r-band (5620 < A < 6980 A) to map out [O ll]A3727A. The last source, 
3C67, is at a lower redshift, so we observed it only in the r-band for the [O ill] coverage. We 
used the R400-G5305 grating leading to a dispersion of 0.687 A per pixel during the z-band 
observations and 0.680 A for the r-band. 

For each object in each band, we took two dithered exposures of 1800 sec each with 
an 0'.'5 offset along the direction of the longer (7") side of the rectangular field of view, 
with the exception of SDSS Jl 101+4004 where only one 1800 sec exposure was obtained. 
The two exposures that go into each observation are reduced separately. We perform the 
data reduction using the Gemini package for IRAF0, following the standard procedure for 
CMOS IFU de^crib^d in the ta=b .n^osinfoifu and gn^osexan^^iefl except that (a) we use 



^The Image Reduction and Analysis Facility (IRAF) is distributed by the National Optical Astronomy 
Observatories which is operated by the Association of Universities for Research in Astronomy, Inc. under 
cooperative agreement with the National Science Foundation. 

^http://www. gemini.edu/sciops/data/IRAFdoc/gmosinfoifu. html 
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an overscan instead of a bias image throughout the data reduction, and adjust the relevant 
parameters so that the bias correction is apphed only once on each image, and (b) we set 
the parameter "weights" of gfreduce to "none" (in contrast to "variance" as suggested by 
the standard example) to avoid significantly increased noise in some parts of the extracted 
spectrum. 

For each dither position, the routine gfreduce trims, overscan-subtracts and extracts 
Gemini Facility Calibration Unit (GCAL) and twilight fiat observations, and gfresponse 
makes response curves with twilight correction. We similarly process the copper-argon arc 
calibration images using gfreduce and then perform wavelength calibration using gswave- 
length and gftransform. We use gfreduce to trim and overscan-subtract science frames, but 
then we interrupt the standard pipeli ne to remove the c osmic rays using the spectroscopy 



version of the L. A. Cosmic software (Ivan Dokkunj|200ll ). The cleaned output is then sent 



back to the master task gfreduce for further processing: the data from separate CCDs are mo- 
saiced and fiat-fielded, the traces of the GCAL flat are referenced for extracting the spectra, 
the wavelengths of the spectra are calibrated, and the derived sky spectrum is subtracted. 
As the final step, the reduced data from each exposure are resampled and interpolated onto a 
data cube with a spatial sampling scale of O'.'l, and the two frames are shifted and combined 
to produce the final science data cube. 



2.3. Calibration 



We flux-calibrate our data using the spectra of our science targets from the SDSS Data 
Release 7§|. SDSS spectra are collected by fibers with a 3" diameter at a typical seeing of 
~2", and SDSS spectrophotometry is better than 10% over its entire wavelength coverage 



of 3900-9100A flAbazaiian et all 120091 ). Since SDSS fiber fluxes are calibrated using point 



spread function (PSF) magnitudes, SDSS spectrophotometry is corrected for fiber losses. 

We subtract the host galaxy continuum from both the SDSS and the IFU observations 
and concentrate on calibrating the [O ill] observations. In order to do this, we first locate 
the pixel (the 0.1" "spaxel") in which the [O ill] line has the broadest wing on the spectrum. 
Then we manually select a pair of wavelength intervals on both sides of [O ill] where the 
continuum is free of any line emission and artifacts (e.g., chip gaps, strong residuals of sky 
line removal). This pair of wavelength intervals is then fixed for all the pixels, and we linearly 
interpolate between them to define the continuum, which we then subtract the best-fit line 
for each pixel. 



■^http: / /www. sdss.org/dr7 
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The seeing of our IFU observations varied between &!4 and &!7 during the three ob- 
serving nights. To simulate the SDSS fiber observations, we convolve the IFU image at each 
wavelength with a Gaussian kernel whose Full Width at Half Maximum (FWHM) satisfies 
FWHM^ + seeing^ = 2"^ to mimic the SDSS observing conditions and then we extract the 
spectrum using a 3"-diameter circular aperture. We then collapse the spectrum within the 
wavelength intervals used for local continuum fitting. The resultant [O ill] intensity, when 
compared to the fiux measured from the SDSS spectrum whose local continuum under [O ill] 
is subtracted in the same way as we do for the IFU data, provides a fiux calibration factor 
for all emission lines. The detailed velocity structure of the [O ill] emission line is highly 
consistent between the 3" extraction from the IFU data (after spatial smoothing) and the 
SDSS fiber spectrum once the SDSS vacuum wavelengths are properly recalculated in the 
air. 

The calibration of the SDSS data includes a PSF correction to recover the fiux outside 
the fibers assuming a 2" seeing, which needs to be removed for our purpose. In the last 
step of our calibration, we downgrade the image of our standard star to 2" resolution and 
calculate the correction factor for a 3" circular aperture centered on the star, which changes 
the final fiux calibration by 7%. We then take this into account for the final calibration of 
the IFU data against the SDSS spectra. 



3. Measurements of the ionized gas quasar nebulae 

3.1. [O III] maps and surface brightness profiles 

We obtain a fiux-calibrated [O ill] line intensity map by collapsing the uncalibrated 
continuum-free IFU datacube over the [O ill] wavelength range and multiplying it by the 
calibration factor described in Section 12. 3[ The [O ill] maps of the eleven radio-quiet type 2 
quasars are shown in Figure [H where the false color is used to represent the intensity on a 
logarithmic scale. We zoom in and out on the objects depending on their redshifts to display 
all images on the same linear scale. We determine the noise level in each pixel on the map 
using wavelengths outside but close to the [O ill] line. The surface brightness sensitivity 
(r.m.s. noise) of our [O III] maps is in the range a = 1.0-2.2 x 10"^'' erg s~^ cm~^ arcsec"^. 
We use a S-cr threshold to create these maps. 

The average seeing of 0.4-0.7" corresponds to a typical linear resolution of ~3 kpc at 
the typical redshifts of our sample. In Figure [2] we show a detailed comparison of the [O 
III] surface brightness profiles with the PSFs. We fit elliptical isophotes to the continuum- 
free [O III] line intensity map of each object using the ellipse task from the IRAF STSDAS 
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SDSS J01 49-0048 



SDSS J0210-1001 



SDSS J0319-0019 



SDSS J031 9-0058 



o O 



SDSS J032U0016 



SDSS J085844417 



3C67 ; J02244-2750 



o 



SDSS J0759t1 339 



SDSS J084U2042 



SDSS J0842+3625 





SDSS J1039'f4512 



SDSS J1040+4745 



O o O 



SDSS J0807+4946 



SDSS J1101 +4004 



\ 



Fig. 1. — Intensity maps of the [O III] line from the sample obscured quasars (11 radio-quiet 
and 3 radio-loud/radio- intermediate objects, separated by the dashed line), with a cutoff at 
a signal-to-noise ratio of 5, shown on a logarithmic color scale. The tick marks of the color 
bar are in units of 10~^^ erg s~^ cm~^ arcsec"^. All maps have the same spatial scale as 
shown in the top le ft panel. The orienta tion of the scattering axis for SDSS J0842+3625 
(101° east of north. IZakamska et al.ll2005l ) is marked by thick dashed lines. The orientation 
of the jets from the radio objects is denoted by thick solid lines. 
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package. In the Figure, we show the surface brightness of each isophote as a function of 
both semi-major and semi-minor axes. The seeing is determined by averaging the directly 
measured FWHM of a sample of field stars in the acquisition image taken right before the 
science exposure. In Figure [21 each PSF is represented by the radial profile of our standard 
star whose width is scaled to match the seeing of each observation. The standard star is 
slightly elongated (ellipticity = 0.07), and here we conservatively use the profile along its 
major-axis. The radial profiles of our targets are clearly extended except SDSS J0841-I-2042 
and SDSS J1039+4512. These two objects, although they appear only marginally resolved, 
demonstrate evident velocity gradient and [O ill] line width variation across their maps 
(which we will show in an upcoming paper) and a strong radial change of [O lll]/H^ line 
ratio similar to that seen in other objects (see Section . We thus conclude that all the 
quasar nebulae are spatially resolved by our IFU observations. 



3.2. Size Measurements 

The physical size of a quasar's narrow-line region is one of the basic parameters that 
characterize the impact of the active nucleus on its environment; yet this value is not an 
easy one to define. The options include (but are not limited to): 

1. i?5cr, the radius of the best-fit ellipse which encloses pixels with a signal-to-noise ratio of 
5 or higher. This measurement faithfully reports the observed extent of a quasar nebula, 
but the size depends on the depth that an observation reaches, making it difficult to 
compare objects observed with different sensitivity and/or at different redshifts. 

2. Res, the effective radius which encloses half of the total luminosity of the line emission. 
Due to the steep decline of the light profiles in the central parts, this approach has 
the advantage that Res is relatively insensitive to the detection threshold or to faint 
extended emission. However, this is also the disadvantage if faint extended emission 
is what we are after. Furthermore, since there is no universal surface brightness pro- 
file amongst quasar nebulae, i?eff does not provide a complete characterization of the 
nebula. 

3. -Robs, the isophotal radius at an observed limiting surface brightness of 10~^^ erg s~^ 
cm~^ arcsec"^. In this work, this threshold is chosen so that the measurements stay 
above ~ 3a sensitivity level for our sample objects. This definition may be useful for 
nearby objects in which cosmological dimming effects are negligible. 

4. -Rint, the isophotal radius at an intrinsic limiting surface brightness of 10~^^ erg s~^ 
cm~^ arcsec"^, which is corrected by a factor of (1 + z)^ (ranging from 4.5 to 7.4 
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Fig. 2. — [O III] surface brightness profiles of the quasars in our sample along both the semi- 
major axes of the best-fit isophotal ellipses (solid lines) and the semi-minor axes (dot-dash 
lines), compared to the profiles of continuum emission along the semi-major axes (dashed 
lines) and the PSF at the time of the observation (dotted fines). The peak [O ill] surface 
brightness (/q) marked on each panel is in 10~^^ erg s~^ cm~^ arcsec"^. The continuum 
emission profiles are truncated at their respective 5-a levels. Also given in each panel is 
the absolute value of the best-fit exponent (77) of power-law fits to the outer regions of the 
nebulae along the semi-major axes (i.e., Ir,io m] oa R"^). 
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for our objects) in order to account for the cosmological surface brightness dimming 
effect. This threshold corresponds to a 6-0" detection for the worst case in our sample. 
This relatively conservative threshold is a compromise between our data and the data 
with lower sensitivity and/or at higher redshifts which we take from the literature for 
various comparisons. i?int is better physically motivated than the previous measures, 
in that it does not depend on the depth of the observations and the redshifts of the 
objects. 

All the above radii are defined as the semi-major axes of the best-fit ellipses at their 
respective isophotes or signal-to-noise {S/N) levels. The measured diameters of the ionized 
gas nebulae range from 15 to 40 kpc, as seen from the R^^^ values determined from the [O ill] 
maps. Not just the faint extended gas envelopes, but the bright central parts are resolved 
by our observations: the effective radii are from 1.4 to 2.9 times the effective radii of the 
corresponding PSFs. All size measurements are reported in Table [2l 

We also compute the ellipticities of the isophotes corresponding to the radii listed above. 
Ellipticities are defined as e = 1 — b/a where a and b are the semi-major and semi-minor axis, 
respectively, and are reported in Table |2] except for ellipticities at i?efr which are affected by 
the PSF and are thus not reported. The nebulae of radio-quiet quasars are nearly round, 
with eso- = 0.13 ± 0.12 (mean and standard deviation). 

Furthermore, we create continuum maps for all objects by summing up the spectrum 
over a 2OO-3OOA wavelength range where line features are absent and artifacts are minimal 
(typically either 4400-4700A or 5200-5400A), and perform differential photometry on these 
maps in the identical manner to that employed for [O ill] line images. We list both the Rl'^^ 
and the R^^^ in Table [2J Our sensitivity to the continuum observations is limited, since the 
continuum is extracted using a relatively narrow wavelength range; therefore, the continuum 
profiles shown in Figure |2] with dashed lines are truncated at much smaller distances from 
the center than the profiles of the nebular emission. The continua are spatially resolved, 
with -Rgg'^* ranging from 1.3 to 2.3 times the effective radii of the corresponding PSFs, and 
are nearly round for the radio-quiet quasars, with el^^ = 0.14 ± 0.05. 

In Figure [3] we show the distribution of relative sizes and relative ellipticities for the 
nebulae and the continua. For radio-quiet objects, the ratio Rcs/Rlfp^ ranges between 0.90 
and 1.33, with the mean and standard deviation of 1.11 ± 0.18. The continuum is likely 
due to the combination of the starlight from the host galaxy and the small amount of 
light from the obscured quasar that escapes along unobscured directions, scatters off the 
interstellar medium of the galaxy and reaches the observer. Since scattered light is polarized, 
polarimetry or spectrop olarimetry can provide a way to determine t he dominant direction 



of quasar illumination (lAntonucci fc Millerl Il985l : iTran et al.l Il995l : IZakamska et al.l 12005 
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20061 ). The contribution of the scattered hght to the total continuum emission ranges from 
to 60% (ILiu et aLll2009l ). Thus the spatial profile of the continuum is determined by both the 
distribution of stars in the galaxy and the distribution of the interstellar medium, whereas 
the sizes of the emission line nebulae are determined by the distribution of the ionized gas. 

The small ellipticities of the continuum emission suggest that th e scattered light emis- 
sion (which tends to show pronounced conical or biconical structure; IZakamska et al.ll2006[ ) 
is either not an important contributor or is concentrated close to the center of the galaxy 
on scales smaller than the seeing. Of the sources presented here, spectropolarimetric obser- 
vations are available for SDSS J0842-I-3625. This object's continuum emission is polarized 
at an astonishing 20% level, and there are no discernibl e stellar features in the deep optical 
spectrum we took of this source (IZakamska et al.ll2005l ). These data indicate that much, if 
not most, of the continuum is due to the scattered light. The position angle of the E- vector 
of polarization is 11° east of north and therefore the dominant illumination direction is 101° 
east of north, as indicated in Figure [1] with dotted lines. Although some elongation of [O ill] 
emission is vaguely visible in this direction, the continuum morphology shows no features. 
Therefore, even in this rather extreme source the scattered light emission is not spatially 
resolved by the ground-based observations, and the combination of HST and polarimet- 
ric observations would be necessary to determine the relative importance of starlight and 
scattered light in our sample. 

Since the narrow-line emission and the continuum are produced by different components 
of the galaxy via different emission mechanisms, the close similarity of the effective sizes of 
the nebulae and the continua is rather surprising. The key question is whether the warm 
ionized gas observed via its line emission is co-spatial with stars or whether the gas is in fact 
being pushed out of the galaxy by the quasar into the intergalactic medium. Therefore, the 
proper comparison would be not between the effective radii of the nebular and continuum 
emission, but rather between the density profiles of the ionized gas and the stars. We hope 
to be able to obtain the density profiles of the ionized gas by extending the range of our 
photo-ionization models of the kind presented in Section 14.41 To get accurate stellar light 
profiles would require _&ST-quality data to separate out the scattered light contribution 
JZakamska et ahlbood . [2006h . 



3.3. Size-Luminosity Relation 

We now investigate the relationship between the integrated [O ill] luminosity and the 
size of the nebula - the size /luminosity relationship shown in Figure |H In addition to our own 
data on the high-luminosity end of the diagram, we use several samples from the literature. 
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These datasets include IFU maps of two ty pe 2 quasars ( Humphrey et al.ll2010l ). long-sht 
spectroscopy of type 2 quasa rs at 2^=0 .1-0. 



of nearby Seyfert 2 galaxies (Bennert et al. 



imaging campaigns ( ISchmitt et al. 



2003 



(jGreene et al. 



2006 



201 11) . and l ong-slit spectroscopy 



Fraquelli et al.ll2003l ). HST narrow-band 



Bennert et al.ll2002l ) are included in the figure (rep- 



resented by gray symbols) but are not used in any quantitative measurements, because they 
are shallower than IFU and long-slit spectroscopy by 1-2 orders of m agnitude and thus d o 
not reach down to the surface brightness levels we use in our analysis (jGreene et al.l 120111 ). 



Except for the solid gray symbols, the radius reported in Figure [His Ri^t, the isophotal 
radius at 10-^7(1 + 

erg s ^ cm ^ arcsec ^ which corrects the sizes for cosmological 
dimming and is thus the measure of choice for samples that include objects at different 
redshifts. For our data, we employ a conservative approach to estimate the uncertainty 
of -Rint- Taking into account that the SDSS calibration error is < 10% (Section 12. 3p . and 
that the isophote at -Rint corresponds to a signal-to-noise ratio of 3.5-10.7 among the radio- 
quiet objects, we find a 10-30% uncertainty. This is an overestimate, because measuring an 
isophotal radius involves averaging many pixels and the uncertainty is thus suppressed; but 
the pixels are actually not independent and the uncertainty is difficult to further quantify. We 
budget a further 30% uncertainty to account for continuum subtraction and instrumental 
and pipeline errors, and combining all these uncertainties in quadrature leads to a final 
uncertainty of < 40% in our surface brightness measurements. Since the surface brightness 
declines as a steep function of radius, we find that the uncertainty in our determined Rmt 
are about 5-14% for the eleven radio-quiet quasars. 



We remeasure the radii of the targets in the four relevant s tudies (IFraquelli et al.ll2003l : 
Bennert et al.l l2006t [Humphrey et al.l l2010l : iGreene et al.l 1201 ll ) to determine Rint in a uni- 
for m manner to construct the size-luminosity relation. Among the six objects observed 
by Humphrey et al.l (120101 ) . only two radio-quiet quasars are resolved and found to be sur- 
rounded by ionized gas nebulae. We assign a 50% uncertainty for their radii, which is an 
arbitrary but conservative estimate fro m their IFU [ O ill] maps with 1" spaxels and 1.2"- 
-1.6" seeing. All the 15 quasars in the iGreene et al.l ( 120111 ) sample are measured using the 
1- dimensional spatial p rofile of [O ill], and we adopt the uncertainty of about 50% estimated 
by IGreene et al.l ( 120111 ) to account for the unknown elongation of these nebulae. 



For the objects in the iBennert et al.l ( 120061 ) and IFraquelli et al.l (12003! ) samples the 
cosmological dimming effect of surface brightness is not important because these Seyfert 2 
galaxies are nearby (the most distant object CGCG 420-015 is at 2; = 0.0294). Our limiting 
surface brightness level corresponds to lO'^^ '' erg s^^ pc^^ in their units for all these targets. 
Both groups have provided the best-fit power laws of the [O ill] surface brightness profiles, 
which we use to obtain Ri^t (the [O ill] profile of ESO 362-G008 is not available and we 
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exclude this source from the sample of iBennert et al.ll2006l ). For the two ta rgets included 
i n bot h samples, NGC 1386 and NGC 5643, we choose the profiles given in iBennert et al. 
( 120061 ). because these authors subtract the contributi on of the star fo r matio n to the ionization 
of [O III]. We ad opt the uncert a inties estimated by iFraquelli et al.l ( 120031 ) for their objects, 
while for the five IBennert et al.l ( l2006l ) objects we propagate the uncertainties in the best-fit 
power indices to Rint- 

The uniformly defined sizes from previous investigations and the addition of our high- 
quality size measurements at high luminosity are combined into a single size-luminosity 
relationship in Figure HI The i?int sizes are strongly correlated with [O ill] luminosity, and 
we fit the relationship with a power law. We perform a linear regression on the well-measured 
l ogari thmic data points, with upper limits and double continuum sources from iGreene et al. 
( 120111 ) excluded. We use linear regression algorithm linmix_err from the IDL Astronomy 
User's Librarjfl, which applies a Bayesian approach using Mark ov chain Monte Carlo to 
calculate the posterior probabilities as described in iKellyl ( 120071 ). By fitting a Gaussian 
profile to the resultant posterior distribution of the slope and intercept which are both close 
to normal distributions, we take advantage of Bayesian methods to reliably calculate the 
uncertainties. Assuming a 15% uncertainty for all L[o m] measurements, we find a best-fit 
linear relation 



log(i?int/pc) = (0.250 ± 0.018) log(Lro iiii/10^2 s~i) + (3.746 ± 0.028). 



The best-fit parameters are consistent within la with those found by iGreene et al. 



( I2OIII ). In particular, we confirm the shallow slope they find (0.22 ± 0.04 ). We caution 



that t he shallow slope hinges in part on the data points at the faint end from IFraquelli et al. 



( I2OO3I ) , which are based on extrapolations from their power-law fit to t he [O ill] profile rather 
than real measurements (see also discussion in lHusemann et al.ll2012[ l. 



We detect extended [O ill] emission in every case, whereas iHumphrev et al.l (120101) con- 
ducted IFU obseryatiori s of type 2 quasars drawn from the same parent sample ( IZakamska et al 
2OO3I : iReyes et al.ll2008l ) and resolved two out of six objects. The most likely reasons for this 
difference are the smaller luminosity of the objects selected by Humphrey et al.l (120101 ) (which 
implies smaller expected sizes, as per our size-luminosity diagram) and the 2-3 times better 
seeing of our observations. 



"'http://idlastro.gsfc.nasa.gov/ 
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3.4. Emission line ratios 



Among the vari ous optical line ratios, [O lll]/H/3 is one of the bes t diagnostics of the ion- 
izatio n mechanism (IBaldwin et alJIlOSlt IVeilleux &: Osterbrocklll987l : lOsterbrock fc Ferland 
20061 ). In order to create the intensity maps of the H/3 emission line, we first need to evaluate 
the contribution of the host galaxy. Unlike the case of [O ill] where a local linear contin- 
uum fit is sufficient, the measurement of H/3 may be biased by the absorption lines in the 
stellar photospheres of the host galaxy. We take the fiber spectrum of each object from the 
SDSS database and decompose the continuu m into a linear s um of the principle components 
derived from pure absorption-line galaxies (IHao et al.ll2005l ). We subtract the best-fit host 
galaxy continuum from each SDSS spectrum and find that the H/3 fiux changes by no more 
than a few percent, with the exce ption of SPSS J0 319-0019 (Figure [5]). This object shows a 
strong post-starburst continuum (IGoto et al.ll2003l ). and subtraction of this component leads 
to a detection of H/3 emission with a S/N ~ 2.5 at the peak of the line (which is otherwise 
undetected in the total spectrum). We attempted to apply this strategy to the IFU data of 
this target in a spatially resolved manner, but could not reliably recover the extended H/3 
emission, because the baseline determination is hindered by the chip gap in vicinity of the 
H/3 line. Therefore, we do not perform spatially resolved analysis of [O lll]/H/3 in this one 
target. In the remaining objects, we create the intensity maps of the II/3 line emission from 
our sample objects in an identical manner as for the [O ill] line. 

To avoid oversampling within the PSF and to increase the S/N ratio, we rebin the 
intensity maps so that the pixel size matches the full width at Icr of the Gaussian profile 
(i.e., pixel scale = 2a) that represents the corresponding observational PSF. In Figure [6], the 
[O lll]/H/3 line ratio is plotted against the [O ill] line surface brightness (S[o mj) and the 
corresponding isophotal radius (semi-major axis from the elliptical differential photometry 
of the [O lll]A5007 map, see Section [3TT]) . In this spatially resolved pixel-by- pixel analysis, we 
have only considered pixels with S/N ^ 3 in both lines. Although the three radio-loud and 
radio- intermediate quasars are included, their isophotal morphology is not as well represented 
by ellipses as radio-quiet objects, therefore the E[o ni]-to-radius conversion introduces larger 
uncertainties. 

Figure [6] reveals a universal behavior of the [O lll]/H/3 line ratio in our entire sample 
of radio-quiet quasars: the ratio persists at a constant value of roughly 10 in the central 
regions, until it reaches a "break" isophotal radius ranging from ~4.1 to 11 kpc where it 
starts to decrease. In Section 14.41 we further discuss this figure and present our qualitative 
interpretation of this behavior. 



In addition to H/3 and [O ill] we detect several weaker emission lines. Because of their 
lower total fiuxes, they are typically detected out to smaller distances than H/3 (and of course 
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smaller yet than [O ill]). One of these, He ll A4686A, turns out to be diagnostic of physical 
conditions in the nebula. We describe the underlying physics and present the He ll A4686A 
measurements in Section! 



The wavelength coverage of the CMOS spectra is rather limited, and H7 is covered 
only in a couple of objects in the region of poor filter transmis sion. The typical reddening 



of the narrow-line regions of type 2 quasars was estimated by (IReyes et al.ll2008l ) from the 
spatially integrated SDSS spectra to be Ay = 1.0 — 2.6 mag based on the observed Ha/H/3 
and H/3/H7 ratios. The only two line ratios we are discussing in the present paper are [O 
lll]/H/3 and He ll/H/3. Because all these lines are close together in wavelength, de- reddening 
would result in a small correction to these ratios of < 0.06 dex. 
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Fig. 3. — Comparison of the extents and ellipticities between the [O IIl]A5007 Une and 
the continuum maps of our quasars. In the eleven radio-quiet objects (striped bars), the 
effective radii of the nebulae and the continua are similar, as are their ellipticities. The three 
radio-loud objects are represented by gray solid bars. 
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Fig. 4. — Size-luminosity relation of the radio-quiet quasars in the [O ill] emission line 
from our sample (black sol i d dia monds), added to the previo us observations incl uding IFU 
mapping ( Humphrey et al. 2010 ) and long-slit spectroscopy ( Greene et al. 201lh of type 2 
quasars at z=0.1-0.4, and l ong-slit spectroscopy of nearby Seyfert 2 gal axies (IBennert et al. 



200d: IPraaudl 



Bennert et al 



et al.l l2003l ). HST narrow-band imaging campaigns (ISchmitt et al 



2003 



200 



^ are shown with gray symbols but are not used in any quantitative 



analyses, because th ey are shallower th an IFU and long-slit spectroscopy by about 1 to 2 
orders of magnitude (iGreene et al.ll201ll ). 




Fig. 5. — Spectra of our sample quasars, taken from the SDSS DR7 except for 3C67 whose 
spectrum is obtained by collapsing our Gemini IFU datacube along spatial dimensions. The 
original SDSS spectra have much broader wavelength coverage than the GMOS spectra 
examined in this paper, but do not contain any spatial information. The marked emission 
lines are: (a) [O II] 3727A; (6) [Ne ill] 3869A; (c) [Ne ill] 387lA; (d) H5 4102A; (e) H7 
4341A; (/) [O III] 4363A; (g) He II 4686A; (h) H/3 486lA; (i) [O III] 4959A; (j) [O III] 5007A. 
Right panels show the velocity substructures in the [O III] 5007A line in the rest frame of 
the sources as determine using the redshifts listed in Table [H The GMOS spectra analyzed 
in this paper typically cover He ll, H/3 and [O III]. 
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Fig. 6. — [O lll]/H/5 as a function of [O ill] surface brightness and isophotal radius (semi- 
major axis). For each object, we hst the [O IIl]/H/3 ratio on the plateau (labeled "const"), 
the [O III] surface brightness and radius where the break occurs, and the best-fit power-law 
slope of the non-plateau regimes for each object are marked. The breaks are depicted by 
dotted lines. 
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Discussion 



4.1. Comparison with radio-loud objects 



We find extended (> 20 kpc-scale) ionized gas nebula in every highly luminous radio- 
quiet obscured quasar that we have ex amined at z ^ 0.5. T h is find ing is in contrast to 
many previous studies, beginning with IStockton &: MacKentyl (119871 ). which found few if 
any radio-quiet objects with extended line emission; rather, it was the radio-loud quasars 
that routinely displayed powerful extended emission-line regions. Since we are now able to 
extend previous studies to a higher quasar luminosity, particularly for obscured quasars, we 
revisit the comparison between the properties of the nebulae of radio-loud and radio-quiet 
objects. 

Jets are frequent ly implicated in strong interactions with galactic int erstellar matter 
(lOosterloo et al.ll2000l ). in exciting the gas through shocks (ICecil et al.ll2001[) and in produc- 



ing outflows of gas from the host gala.xy, of t en on scales of severa l kpc ( iBaum et al.lll990 



19921 : iMcCarthvl Il993l : iMoreanti et all l2005l : iNesvadba et al.ll2006f ). At high redshifts, the 
emission-line regions can be extended over tens of kpc, contain significant amounts of mas s 
and are typically well-aligned with the large-scale jet emission (INesvadba et al.ll2006l . |2008| ). 
At low r edshifts, there is no particularly strong alignment between the jet and the emission- 
line gas (jPrivon et al.ll2008l ). perhaps with marginal evidence for alignment in younger radio 
sources still confined to their host galaxy. The lack of alignment could have many explana- 
tions - for example, the pattern of quasar illumination of the interstellar medium may be 
different from the jet direction, or inhomogeneities in the gas could lead to nebulae expand- 
ing in the directions of lowest density. No obvious alignment is seen in the three radio-loud 
sources we observed as part of our campaign (bottom of Fig. [1]). Of the four extended 
emission features seen in our three radio-loud sources, only one feature in SDSS J0807-I-4946 
lies along its radio axis. 

For comparison with our objects, an ideal radio-loud sample would contain type 2 sources 
(i.e., radio galaxies) with similar redshifts and [O ill] lumin osities observed to sirn ilar depth. 
Narrow-band imaging data from space (e.g., for 3C objects. iTremblay et al.ll2009l ) are there- 
fore not a good option because they are shallower t han our IFU data by 1-2 orders of 
magnitude. A close-to-optimal sample was studied by iFu fc Stockton! (120091 ) who consider 
eight radio-loud sources. Their objects are type 1 radio-loud quasars which are a bit closer 
than ours (z = 0.16 — 0.37) and have slightly lower luminosities (L[o m 
but their IFU data is very well matched to our observations. 



1042-43 g-l^ 



The r nost striking difference between our radio-quiet objects and the radio-loud sources 
studied by lFu fc Stockton! (120091 ) is that the objects from the latter sample show much more 
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complex morphologies in their [O ill] emission. The nebulae studied by lFu &: Stockton! (j2009[ ) 
often consist of several kinematically and/or morphologically separate clumps and blobs, and 
even when the emission is connected into a single nebula it is not a linear structure that we 
could fit with an ellipse for comparison with our observations. iFu fc Stockton! ( 120091 ) suggest 
that the blobs and clumps could be due to quasar-illuminated merger debris or the remnants 
of a quasar-driven wind. 

In view of the limited availability of comparison radio-loud objects at low redshifts, we 
consider high-redshift samples as well. We construct the radio-loud comparison sample by 
combining the three radio-loud quasars in our own Gemini campaign (Fig. (H bottom three) 
with three powerful radio galaxies MRC 0316-257, MRC 0406-244 and TXS 0828+193. These 
thr ee objects, locat e d at z ~2-3, were observed with SINFONI in the IFU mode on the VLT 
by iNesvadba et al.l (12008! ). The isophotal radii i?int (semi- major axes) and ellipticities eint 
are measured the same way as for radio-quiet quasars at the [O ill] surface brightness of 
]^g-i5 gj,g g-i Qj-Q-2 arcsec"^ (i.e., corrected for cosmological dimming). These results, along 
with some basic information about the three high-redshift objects, are listed in Table |3l 

In Figure [HI we show the radii -Rint and ellipticities eint of the eleven radio-quiet and 
the six radio-loud quasars plotted against their radio luminosity at 1.4 GHz. The 1.4 GHz 
fluxes a re taken from the FI RST survey for SDSS objects and the NRAO/VLA Sky Survey 
(NVSS, ICondon et al.l Il998! ) for the rest and are then used to ca l culat e the rest-frame- 
corrected monochromatic luminosities at 1.4 GHz (jZakamska et al.! 120041 ). In the case of 
non- detections in the FIR ST survey, we report the upper limits for a point source (5(7+0.25 
m Jy, ICoUinge et aP l2005h . 

Separated by the vertical dotted line at z/Ly[1.4GHz] = 1.1 x 10"^^ erg s~^ (our adopted 
radio- quiet /radio- loud criterion discussed in Section [^7T]) . the nebulae of the two populations 
of quasars demonstrate an evident difference in their ellipticities. All the radio-loud objects 
have elongated narrow-line nebulae with eint > 0.25, while ten out of eleven radio-quiet 
quasars have eint < 0.25, with SDSS J0321+0016 being the only exception. In terms of their 
respective mean value and standard deviation, the ellipticities of radio-quiet quasars have 
Cint = 0.16 ± 0.10, while their radio-loud peers have e^nt = 0.45 ± 0.14. The most striking 
morphological difference between the ionized gas nebulae around radio-loud and radio-quiet 
objects is that the radio-quiet ones are significantly round er. In contrast, the [O ill] line 
emission from all the three objects in lNesvadba et al.! (120081 ) are remarkably elongated along 
the direction of their radio jets. 

The semi-major axes of radio-loud nebulae (-Rint = 19-9 ± 6.9 kpc) are larger than those 
of radio-quiet objects (11.1 ± 3.0 kpc), although the difference in sizes is not as dramatic 
as it is for ellipticities and there is some overlap in the size distributions. It is perhaps the 
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larger sizes and larger ellipticities of the nebulae in radio-loud objects that made them more 
easily detectable in previous studies. 

Our radio-loud / radio-quiet criterion is set at vL,j[lA GHz] = 1.1 x lO'*^ erg s~^, which 
is a negligible fraction of the bolometric luminosity, but could cor respond to a significant 



mechanical energy of the jet. For example, applying equations of ICavagnolo et al.l (120101 ) 
we find that an Lmoch = 5 x 10^^ erg s~^ jet could be hiding behind this rather small radio 
luminosity. It is precisely the dramatic difference in morphology and ellipticity of the ionized 
gas nebulae that occurs at roughly this radio luminosity that leads us to conclude that jets are 
not an important contributor to supplying and exciting ionized gas in our eleven radio-quiet 
sources. 



4.2. Comparison with Unobscured Quasars 



Another natural comparison is between our sample and unobscured radi o-quiet quasars 
atsimilar lu minosities and redshift s. Among the existing studies, the sample of lHusemann et al, 
( 120081 ) and (iHusemann et al.ll2012[ ) is fairly well matched to ours in luminosity, although their 
sources are at slightly lower redshifts [z < 0.3). At similar surface brightness limits, those 
authors do not find the spatially extended round nebulae that we see in our radio-quiet 
sample of type 2 quasars. They find extended [O ill] emission in 61% of the objects in 
their sample (the rest are non-detections). Compared to our objects, their nebulae appear 
elongated and/or irregular (with ellipticities range from 0.15 to 0.64), and there is a trend 
for targets with high Fe ll equivalent width (the doublet AA4924,5018 plus the complex at 
5100-5405A) and low H/3 FWHM to have smaller nebulae. This last issue is not testable 
with our data because the Fe ll emitting regions are obscured in our targets. These authors 
also find that the extended nebulae are aligned with the radio jets even when the jets are 
weak, and that the size of the nebulae increases for increasing radio luminosity. 

One of the possible expla nations for the di f ferenc e between the detection rate of nebulae 
in our sample and in that by iHusemann et al.l (|2008[ ) is that by selecting sources with high 
[O III] equivalent width (as is the case in our sarnple) we are preferentially selecting those 
with high covering factor in O^"*" ( iBaskin Sz Laoiil2005t iLudwig et al.l 120091 ). Furthermore, 
despite the careful subtraction of the bright PSF (the quasar itself), faint emission on the 
scales of the host galaxy may be difficult to see in type 1 quasars, and thus the detected 
nebulae may be biased toward those that are larger and more elongated. Finally, geometric 
effects of quasar illumination may play a role: since type 2 quasars are obscured along the 
line of sight, the illumination of the gas may be happening on average closer to the plane of 
the sky than it does in type 1 quasars, thus producing apparently larger nebulae. 
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Alternatively, the difference in the morphology of ionized gas around type 1 and type 
2 quasars may reflect an actual evolutionary difference between the two types. Theoretical 
models of galaxy formation have long postulated an evolutionary scenario in which galaxy 
mergers induce both star formation and nuclear activity, triggering a transition from an 
obscured accretion and star formation stage (more likely to be s een as a type 2 source) to an 
unobscured phase as a type 1 quasar (e.g.. ISanders et al.lll988l ). A variety of observational 
studies support a scenario of this type, finding for example that host gal axies of type 2 active 



galac t ic nuclei exhibit more star formation than do the hosts of type I's ( lHdl2005t iLacy et al. 



20071 : IZakamska et al.ll2008l ). Thus it is plausible that obscured quasars are more likely to 



be seen with extended ionized emission because they are in fact asso ciated with a phase in 
which gas is being expelled out from the galaxy jHopk ins et al.ll2006l ). 



4.3. Comparison with Seyfert galaxies and the origin of gas 



Schmitt et al.l (120031 ) compiled and analyzed HST narrow-band imaging data on [O ill] 
nebulae in a large sample of Seyfert 1 and Seyfert 2 galaxies. The typical r.m.s. surface 
brightness limit of their observations is 1.2 x 10~^^ erg s~^ cm~^ arcsec"^, so these data 
probe the bright central regions of the targets on scales of a few hundred parsecs at typical 
resolutions of ~50 pc. The median [O ill] luminosity of the Seyfert galaxies in this sample 
is 2.5 dex smaller than the median [O ill] luminosities of the quasars in our sample. 

Taking just the 38 type 2 Seyfert galaxies from this sample, we find that the ellipticities 
of their nebulae are distributed between 0.06 and 0.75, with the mean of 0.42. Many objects 
show pronounced biconical structures in their [O ill] emission lines. Fu rthermore, we f ind 10 



objects with position angles of the [O ill] nebulae ^[o m] measured by ISchmitt et al.l (120031 ) 



for which pola.rizati on position angles 6'poi are available in the literature (as referenced in 



Schmitt et al.l 120031 ). The average difference between the two sets of angles is (|A6'|) = 
97° ±6°. This value is in strong disagreement with that expected if 6\o m] and 6'poi were each 
independently drawn from a uniform distribution between 0° and 180°, in which case the 
probability distribution of |A6'| would be p{A9) oc (1 — |A6'|/180°) and would be peaked at 
small angles. 



As was already discussed in Section [3l2| the polarization is produced when light from the 
hidden nucleus is scattered off of the extended material, and the illumination direction (as 
projected on the plane of the sky) is orthogonal to the measured polarization angle. Thus, 
in Seyfert 2 galaxies [O ill] nebulae are not only elongated, but are closely aligned with 
the illumination direction. Both the scattered light and the [O ill] emission are produced 
along the directions of clear view toward the nucleus, and in the classical model of toroidal 
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obscuration (lAntonucci fc Millerl Il985l) one ex pects to see ionization bicones - as indeed 
observed in the sample of ISchmitt et al.l (|2003[ ). 



In our sample, the observed morphologies of the nebulae are determined both by the gas 
distribution in and around the host galaxy and by the illumination pattern of the quasar. The 
lack of elongation or biconical structures either in the morphology of the ionized gas or in the 
continuum tracing the dominant quasar illumination directions is one of the major surprises 
of our observations. One possible explanation is that the typical opening angles of quasar 
illumination are large, producing quasi-isotropic illumination pattern. Another possibility is 
that at the high luminosities explored in this paper the geometry of obscuration is different 
from that seen in Seyfert galaxies: perhaps there is no organized toroidal structure, but 
rather obscuration is quasi-spherical and patchy, so that ionizing emission escapes along 
random beams (but not toward the observer) which we cannot spatially resolve by our data. 
Both these possibilities are at odds with our HST imaging of somewh at less luminous type 2 
quasars which clearly shows ionization cones in at least some objects (jZakamska et al.ll2006l ). 
although a strong luminosity dependence of the geometry of obscuration would reconcile 
these results. HST imaging of the objects in this sample in the blue band sensitive to the 
scattered light emission can settle this issue. 

An alternative possibility is that shock excitation rather than photo-ioniza tion is respon- 



sible for the observed st rong line emission. The radiation pressure of the quasar (jMurrav et al 



1995 



2012 



Proga et al.lboooh driv es a low-density high-temperature wind (iFaucher-Giguere fc QuataertI 



Zubovas fc Kingll2012n which slams into dense clumps in the host galaxy, resulting in 



shocked emission. Numerical simulations demonstrate that even if there are anisotropies in 
the initial driving of the outflow, the outflow isotropizes over the range of scales involved 
between the scale on which the outflow is accelerated (< 10 pc) and the scale on which it 
sweeps up the interstellar medium of the galax y 10 kpc), prod ucing structures qualita- 
tively similar to those seen in our observations (INovak et al.ll201l[). Shocks propag ating at 
several hundred km s~^ normally do not produce [O lll]/H/3 ~ 10 (IRich et al.l 1201 ll ) as seen 
in our quasars, and thus they are normally disfavored as the primary origin of narrow line 
emission. Furthermore, the efficiency of con verting mechanical en ergy into line emission is 
expected to be low at such shock velocities (INesvadba et al.ll2008l ). One possibility is that 
the shock blast propagates much faster (> 1000 km s~^) than was previously considered, 
resulting in higher [O lll]/H/3 ratios consistent with those observed. Thus, we will continue 
to explore this possibility in our future work since photo-ionization does not at the moment 
provide a fully satisfactory explanation of all of the aspects of our data. 



- 28 - 



U 

Q_ 



30 
25 

20 

15 

10 



0.6 - 



0.4 - 



0.2 - 



0.0 - 



<-• 



<-• 



M 



M 



H 



M 



IXI 



1X1 



IXI 



■ M-. 



radio-quiet quasars • 
rodio-loud/intermediate quasors M 
Nesvadba et ol. (2008) ixi 



39 40 41 42 43 44 

log Lj 1.4GHZ / erg s"^ 



45 



Fig. 8. — Sizes and ellipticities (-Rint and ejnt, see Section 13.21 for their definition) of the 
[O III] nebulae around 11 radio-quiet quasars from our observations, compared to those of 
the six radio-loud objects (Table [3]). The comparison sample of six contains three low- 
redshift radio galaxies from our own observations and three high-redshift radio galaxies from 
Nesvadba et~ZI J2008h . 
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4.4. Physical conditions in the nebulae 



In this Section we discuss a plausible model of a nebula under the assumption that 
the gas is illuminated and photo-ionized in a wide-angle manner, either because of the 
large opening angles in the obscuring material or because of the quasi-spherical patchy 
illumination. Unless noted otherwise, the arguments in this section apply both to a static 
nebula and to an outflow. Throughout this section we use r to denote the three-dimensional 
distance from the quasar and R to denote the distance as seen projected onto the plane of 
the sky. 

We aim to determine plausible physical conditions for the line-emitting gas that we 
observe. Regardless of the mechanism by which the gas reached the observed radii, it is 
reasonable to draw an analogy with galaxy-scale winds, in which there is also observed warm 
line-emitting gas at large radii. In these cases, the gas is multi-phase, with different phases 
of the medium observable in different domains of the electromagnetic spectrum. By analogy 



with the winds driv en by supernova explosions in star-burst galaxies (IHeckman et al.lll990 



Veilleux et al.l 120051 ). the simplest model involves two components. One consists of hot, 
low-density, volume-filling potentially X-ray-emitting plasma, whereas the other is warm 
(T ~ 10'^ K), higher density component concentrated in clumps, shells or filaments (hereafter 
"clouds"). It is the latter that produces the emission lines that we see in optical spectra. If 
the two components are in rough pressure equilibrium with each other, their densities and 
temperatures are related via nhot^l 



hot 



n 



T 

warm-' warm- 



The [O lll]/II/3 ratio shows a remarkably uniform behavior as a function of distance from 
the center across our sample. As Figure [H] demonstrates, for all the radio and radio-quiet 
objects that we observe (except for SDSS J0319-0019 where we have only a weak detection 
of H/9), the ratio remains constant at [O lll]/H/3 = 12.3 ± 2.7 out to distances of 7.0 ± 1.8 
kpc, and declines as log [O lll]/H/3 oc i?-3-3±i.i larger distances. In addition, most objects 
show a characteristic excess (~0.2 dex) of [O lll]/H/3 right at about the break radius. 

We consider two physically distinct models that can qualitatively account for the be- 
havior of the line ratios. The first model capitalizes on the fact that the [O lll]/H/3 ratio 
is strongly dependent on the ionization parameter U = $/47rcr^nH(r), where $ is the pho- 
ton luminosity of the qua sar above the ionization energy of hydrogen of 13.6 eV. Following 
Villar-Martin et al.ll2008l we consider a simple case of clouds of constant density nn = 100 
cm~^ illuminated by a quasar with bolometric lur ninosity Lhni = 10^^ erg s~^. In this model, 
for a typical quasar spectral energy distribution ( iFerland et al.l 119981 ) the ionization param- 
eter as a fu nction of distance is log ?7 = —0.87 — log(r/kpc)^. We then use calculations 
presented in IVillar-Martin et al.ll2008l to determine the [O lll]/H/3 ratios of such clouds (Fig- 
ure IHl left panel labeled "Model 1"). As long as the clouds remain "ionization-bounded" 
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at all distances from the quasar (i.e., the number of 0^^-ionizing photons is insufficient to 
ionize the entire cloud) and the value of U stays above logf/ ~ —2.3, the [O lll]/H/3 ratio is 
constant and close to 10. As r increases and U declines below logf/ ~ —2.3, the [O lll]/H/3 
ratio declines as well, in qualitative agreement with our observations. 

Despite the apparent success of this model in reproducing the [O lll]/H/3 ratio, the 
input physics is not entirely satisfactory. In particular, since the clouds' thermal balance is 
maintained by photo-ionization and recombination, their temperatures are nearly constant 
under a wide range of conditions with T = f x 10'^ K where the factor / is typically between 



1 and 2 ( lKroliklll999[ ). Since the model assumes a constant density inside the clouds, this 
implies that the cloud pressure is independent of the distance from the quasar as well. If 
the clouds are in pressure equilibrium with the confining low-density medium, this would 
imply that the pressure of the entire quasar nebula is constant as a function of radius. This 
is hard to reconcile with models, as well as with observations of starburst-driven winds in 



which pressure declines outwards (IHeckman et al.lll990l ). 



Therefore, we consider another simple model in which the pressure in the nebula follows 
P(r) oc r~^. The clouds are in confined by the pressure of the hot low-density component, 
and since their temperature is almost constant, the density inside them drops as nH(r) oc 
as well, which implies that the ionization paramet er is independent of the distance from the 



quasar. We use the photo-ionization code Cloudy ( iFerland et al.lll998l ) to model emission of 
clouds with a hydrogen density of = 100 cm~^ at r = 1 kpc from a Lboi = lO'*^ erg s~^ 
quasar. These parameters correspond to an ionization parameter of logf/ = —0.87. 

In Figure m center, we show the [O lll]/H/3 ratios of such clouds as a function of their 
physical depth. Large clouds are ionization-bounded (i.e., have a large optical depth) and 
show the familiar ratio of [O lll]/H/3 ~ 10. Smaller clouds, however, allow for the photons 
with energies above 54.9 eV (the energy required to ionize O^^ into O^"*") to penetrate through 
the entire cloud, destroying most O^"*" ions and resulting in much lower [O lll]/H/3 ratios. 
If a cloud is participating in the outfiow, it expands (linear size S oc r^/^), its density falls 
(oc r~^) and its optical depth declines (oc r~^/^). Therefore, as the clouds propagate out with 
the quasar wind, they are expected to transition from ionization-bounded (partly ionized) 
to matter-bounded (fully ionized). The same transition occurs in a static nebula as long as 
clouds have the same mass distribution at all distances. As a result, the [O lll]/H/3 ratio is 
expected to fall as a function of distance from the quasar, in qualitative agreement with our 
observations. 

All the arguments presented thus far apply equally to a static nebula or a wind, as long 
as the pressure profiles are similar. We now consider specifically an outfiowing wind. In the 
matter-bounded regime the B./3 luminosity of each cloud Ldoud oc n^S^ oc r~^. If the number 
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density of clouds falls off as (as expected in a steady-state constant velocity wind), then 
the H/3 luminosity density falls off as r~^. As a result, the surface brightness of the nebula 
is expected to fall off as R~^, which is close to the values seen in our sample (Table |21 the 
mean value of the exponent is —3.5 ± 1.0). 

The major difference between Models 1 and 2 is in the degree of ionization of clouds at 
large distances from the quasar. In Model 1, the ionization parameter decreases outward and 
the clouds have a smaller degree of ionization than they do in the inner parts. In Model 2, 
the hard radiation can penetrate through the entire cloud, ionizing O^"*" to higher ionization 
states. Coinci dentally, the ionization energy of O^"*" (54.9 eV) is almost identical to the 
ionization energy of He"*" (54.4 eV). Thus the same photons that destroy O^"^ and lead to 
the decrease of the [O lll]/H/3 ratio in Model 2 produce He^"*". This leads to an increase in 
He^"*" recombination emission, including one of the transitions in the recombination cascade 
observed as the optical line He ll A4686A. Thus, Model 2 predicts a much higher ratio of 
He ll/H/3 in the outer parts in the nebula than does Model 1. The model predictions are 
illustrated in Figure [H right. The models do not take into account projection effects - along 
any given line of sight, regions at different distance r from the quasar contribute to the 
emission. Nevertheless, it is clear that high He ll/H/3 ratios, especially in combination with 
low [O lll]/H/3, are characteristic of the high-ionization matter-bounded regions and can only 
be produced in Model 2. 

In Figure (TD] we present the measurements of the He ll/H/3 line ratios in the 10 sources 
where the He ll A4686A transition is covered by our data. The He ll line is weaker than H/3 
and is thus cannot be measured out to the same distances from the quasar, although typically 
we are able to detect He ll at distances comparable to or larger than the break radius. There 
is a general tendency for He ll/H/3 ratio to increase outward, and in 9/10 objects it goes 
above log [He ll/H/3]=— 0.5 - too high to be compatible with Model 1. The key upper 
left region of the He ll/H/3 vs [O lll]/H/3 diagnostic diagrams is particularly observationally 
difficult because it is reached only in the outermost parts of the nebulae where [O ill] is 
faint; nevertheless, the line ratios in two objects - SDSS J0210— 1001 and SDSS J0842+3625 
- extend into this region. 

We will continue to refine Models 1 and 2 as we investigate other aspects of our data. 
From the simple model setup presented above, we conclude that the pressure in the nebulae 
is most likely decreasing outward and that line-emitting clouds transition from ionization- 
bounded to matter-bounded as the distance from the quasar increases. 
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Fig. 9. — Left: In the model where all clouds are ionization-bounded and have the same 
density {nu — 100 cm~^), the [O lll]/H^ ratio shows a behavior qualitatively similar to that 
observed in our sources (plotted here as a function of three-dimensional distance r). Center: 
A model with declining pressure profile uses the fact that the [O IIl]/H/3 ratio changes 
dramatically as the clouds transition from being ionization-bounded to matter-bounded. 
Here, we show the [O IIl]/H/9 ratios of clouds of varying thickness at a distance of 1 kpc 
from a 10^^ erg s^^ quasar calculated with Cloudy. Right: The two models predict similar 
[O lll]/H/3 behavior as a function of distance, but different He ll/H/3 ratios (Model 1 shown 
with dashed hne. Model 2 shown with solid hue) . 
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Fig. 10. — He ii/H^ as a function of [O ill] isophotal radius (left panels) and [O lll]/H/3. 
Large values of He ii/H^ (log He ii/H^ > —0.5) strongly suggest presence of matter-bounded 
clouds, especially in combination with low [O lll]/H^ (log [O lll]/H^ < 0.5), as seen in 
SDSS J0210-1001 and SDSS J0842+3625. 
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Continued. 
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Conclusions 



In this paper we present optical seeing-limited integral field unit observations of ionized 
gas around z ~ 0.5 obscured luminous quasars. We use these data to determine the spatial 
distribution of the [O lll]A5007A and H/3 line emission. We spatially resolve the emission 
line nebulae in every case and find that the [O ill] line emission from gas photo-ionized 
by the hidden quasar is detected out to 14 ± 4 kpc from the center of the galaxy (sample 
mean and standard deviation are given throughout this section). Ionized gas nebulae around 
radio-quiet obscured quasars display regular smooth morphologies, in marked contrast with 
nebulae around radio-loud quasars of similar line luminosities which tend to be significantly 
more elongated and/or lumpy. Surprisingly, no pronounced biconical structures expected in 
a simple quasar illumination model are detected. 

The main measurements we present in this paper are (1) the spatial distribution of the 
emission-line surface brightness, (2) the spatial distribution of the [O lll]/H/3 ratios, and 
(3) the relationship between the size of the nebula and the emission line luminosity. When 
fit to a power law, the surface brightness of H/3 declines as i^-3.5±i.o ^ function of the 
distance from the quasar, whereas [O ill] surface brightness declines as R~'^■^^'^^^'^ . The [O 
lll]/H/3 ratios show a remarkably uniform behavior across the sample: the ratio remains 
constant at 12.3 ±2.7 and at i? > 7 ± 2 kpc declines as roughly a power-law function of 
the distance from the quasar with an index —3.3 ± 1.1. Finally, we supplement our sample 
with 31 objects in the literature, define a redshift- and sensitivity-independent nebula size 
and present a size-luminosity relationship uniformly measured over six orders of magnitude 
in line luminosity. The nebular size slowly increases with luminosity oc 2,o-250±o.oi8_ 

The galaxy-wide narrow-line emission regions around our quasars are large, but the 
origin of this ionized gas is unclear. In principle, the quasar may be illuminating pre viously 
existing gas in the galaxy halo (e.g. from previous star formation, ISteidel et al.ll2010l ): or the 
gas may be debr is associated with recent merger activity a nd currently illuminated by the 
quasar (see, e.g. IFu fc StocktonI l2009t IStockton et al.ll2006l ). The remarkable roundness of 
the nebulae and the lack of dis tinct blobs or c lumps is in direct contrast to the morphologies 
seen in radio-loud objects (c.f. IFu et al.ll201ll ). It is hard to understand how tidal debris or 
gas accretion could leave such a smooth and isotropic distribution of gas. One option is that 
the gas is being expelled from the galaxy, presumably driven by the radiative pressure of 
the quasar itself ( iMurray et al.lll995l : iProga et al.ll2000[ ). We cannot rule out that previous 
generations of star formation have left behind the gas now illuminated by the quasar, but 
we do note that the smooth featureless morphologies of the [O ill] nebulae in our sample a re 
similar to those seen in wide-angle outflows recently observed by lRupke &: Veilleuxl (120111 ). 



We put forward a simple model of a clumpy gaseous nebula, in which the lower density 
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hotter component is largely invisible, but the higher-density warm component is the one we 
observe via line emission. In our fiducial model, the gas pressure inside the nebula and in the 
clouds drops off as and the temperature of the emission-line clouds is nearly constant at 
a few xlO^ K determined by the photo-ionization and emission balance. The clouds closer to 
the quasar are denser and are ionization-bounded (i.e., the photo-ionizing photons from the 
quasar are all absorbed on the quasar-facing side of the cloud and cannot penetrate all the 
way), whereas the clouds beyond r > 7 kpc are lower-density, larger and matter-bounded 
(i.e., the recombination rate cannot keep up with the ionization rate, and these clouds are 
fully ionized) . Unlike the model with a constant pressure profile, the model with a declining 
pressure profile quahtatively explains our [O lll]/H^ and He ii/H^ hne ratio observations and 
the surface brightness distributions, and will be tested and refined further by our on-going 
observations and modeling. 

The kinematic analysis of extended line emission around quasars is on-going and will be 
presented in our upcoming papers. We will demonstrate that the nebulae whose morpholog- 
ical properties are discussed here are likely due to ionized gas in the outflow from the quasar. 
If these outflows are found to carry significant kinetic energy and/or momentum (~ 1% of 
the bolometric output of the quasar), they will be excellent candidates for the long-sought 
feedback phase of quasar evolution. 
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Table 1. Properties of our Gemini quasar sample. 



Object name 
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(10) 




SDSS J014932.53-004803, 


,7 


RQ 


1.03 


0.566 
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1800x2 


0.52 


42.87 


0.13 


5, 


,70±0. 


,14 


SDSS J021047.01-100152, 


,9 


RQ 


<1.20 


0.540 
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0.58 
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.OOibO. 


,02 
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,7 
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1.95 
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,20 
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SDSS J032144.11+001638, 
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Notes. - (1) o bject name. (2) R a dio loudness (RQ: radio quiet; RL: radio l oud; RI: radio interm ediate). (3) Radio flux at 1.4 GHz in mjy, taken from the 
FIRST survey jBecker et "aLlll995l : fwhite et aLlllQQTf) and t he N VSS survey jCondon et al.lll99i . for 3C6 7 only). Upper limit s are r eported for undetected 
sources (see text). (4) Redshift, from lZakamska et alj | |2003|) and lReves et al.l ||2008| ) for SDSS objects and lEracleous fc Halperiil ||2004| ) for 3C67. (5) Position 
angle of the field of view (in degrees), measured by the shorter axis from north to east. (6) Exposure time (in seconds) and number of exposures. (7) Seeing 
at the observing site (FWHM, in arcseconds). ( 8) Total luminosity of the [O lll]A5007A line (erg s~^), derived from our data calibrated against SDSS spectra, 
except for 3C67 whose [O III] flux is tak en from [Celderman &: Whittle! l ll994l ). (9) Relative uncertainty of the [O III] luminosity. The uncertainty for 3C67 is 
taken from [Celderman &: Whittle! ll 19941 ). (10) Absolute value of the best-fit power-law exponent of the outer part of the [O III] profile along the major axes 
(i.e., /fljoill] °<: R~^, see Figure[2ll. 



Table 2. Properties of our Gemini quasar sample (continued). 





pcont 


,cont 






Ijcont 
-"eff 




-"-obs 


^obs 


^int 




const 


-n-br 




slope 






s 
















K') 


{°) 




















[It)) 




SDSS J0149-0048 


6.3 


0.06 


9.1 


0.10 


3.1 


2.9 


7.4 


0.12 


6.9 


0.05 


31.8 


4.1 


-4, 


.79±0, 


.37 


3. 


,27±0, 


,11 


SDSS J0210-1001 


4.8 


0.21 


17.7 


0.07 


3.0 


4.0 


17.8 


0.14 


14.4 


0.12 


10.0 


7.5 


-3 


.09±0, 


.22 


2. 


,08±0. 


.14 


SDSS J0319-0019 


9.9 


0.13 


7.6 


0.25 


3.0 


2.7 


7.9 


0.22 


6.7 


0.23 


















SDSS J0319-0058 


5.5 


0.16 


11.5 


0.04 


3.3 


3.5 


11.2 


0.05 


9.6 


0.14 


11.9 


7.5 


-5, 


.32±1, 


.34 


4. 


,64±0. 


.22 


SDSS J0321+0016 


8.2 


0.17 


18.3 


0.42 


4.3 


5.2 


19.3 


0.38 


16.6 


0.42 


16.2 


11.0 


-2 


.66±0, 


.35 


5. 


,02±0, 


,23 


SDSS J0759+1339 


6.0 


0.22 


14.1 


0.02 


3.0 


3.5 


15.4 


0.07 


12.9 


0.08 


16.1 


7.5 


-2, 


.91±0, 


.31 


3. 


,06±0, 


,08 


SDSS J0841+2042 


5.5 


0.07 


11.9 


0.09 


2.8 


2.5 


11.8 


0.10 


10.0 


0.17 


15.7 


6.4 


-3 


.53±0, 


.28 


3. 


,65±0. 


,08 


SDSS J0842+3625 


6.5 


0.17 


15.1 


0.12 


3.1 


3.2 


14.4 


0.18 


13.1 


0.19 


10.3 


9.0 


-4, 


.14±0, 


.56 


4. 


,82±0. 


,19 


SDSS J0858+4417 


6.9 


0.12 


11.7 


0.08 


3.0 


3.1 


11.5 


0.12 


10.6 


0.10 


11.9 


5.6 


-3 


.88±1, 


.16 


3. 


,24±0. 


,16 


SDSS J1039+4512 


4.9 


0.19 


12.2 


0.24 


2.5 


2.6 


11.4 


0.22 


9.3 


0.12 


12.2 


5.8 


-3 


.90±0, 


.65 


4. 


,54±0. 


,14 


SDSS J1040+4745 


6.9 


0.09 


14.5 


0.04 


3.6 


3.4 


14.5 


0.08 


13.0 


0.18 


8.7 


7.6 


-2 


.45±0, 


.75 


2. 


,88±0, 


,20 


3C67/J0224+2750 


5.2 


0.07 


12.2 


0.42 


2.0 


2.3 


19.4 


0.42 


10.1 


0.26 


4.6 


8.3 


-2 


.38±0, 


.90 


2. 


,99±0, 


,18 


SDSS J0807+4946 


5.8 


0.22 


19.4 


0.44 


2.7 


3.3 


19.2 


0.47 


14.6 


0.45 


13.8 


6.3 


-1 


.45±0, 


.45 


3. 


,60±0. 


,04 


SDSS JllOl+4004 


4.7 


0.05 


18.9 


0.26 


3.4 


5.2 


19.3 


0.31 


18.7 


0.31 


12.3 


8.2 


-2, 


.90±0, 


.36 


2. 


,16±0. 


,03 



Notes. — (1) object name. (2 3) Scini-major axis (in kpc) and ellipticity (= 1 — b/a) of the best-fit ellipse which encloses pixels with S/N ^ 5 in the continuum 
map. (4-5) Semi-major axis (in kpc) and ellipticity of the best-fit ellipse which encloses pixels with S/N ^ 5 in the [O lll]A5007A line map. (6-7) Half-light 
isophotal radius (effective radius) of the continuum and [O lll]A5007A line emission (semi-major axis, in kpc). (8-9) Isophotal radius (semi-major axis, in kpc) 
and ellipticity at the observed limiting surface brightness of 10~^^ erg s~^ cm^^ arcsec"^. (10—11) Isophotal radius (semi-major axis, in kpc) and ellipticity 
at the intrinsic limiting surface brightness (corrected for cosmological dimming) of 10~^^ erg cm~-^ arcsec"'^. (12) Mean value of the plateau part of the 
[O lll]/H/9 radial profile. (13) Break radius (semi-major axis) of the [O lll]/H/3 radial profile. (14) Best-fit power-law slope of the [O ni]/H/3 radial profile at 
R > ilbr- (15) Absolute value of the best-fit power-law exponent of the H/3 radial profile at il > ilbr (i.e., Ir^h^ <x. R~^). 
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Table 3. Properties of the radio galaxies from iNesvadba et al.l ( l2008l ). 



Object name 


Type 


/l.4GHz 


z 


Seeing 


logi[0 III] 




toff 


^int 


Cint 








(mJy) 




(") 


(erg s^i) 


(kpc) 




(kpc) 




MRC 0316' 


-257 


RL 


447 


3.127 


0.7x0.4 


44.29 


4.3 


0.37 


21.27 


0.59 


MRC 0406- 


-244 


RL 


626 


2.425 


0.7x0.6 


44.50 


4.6 


0.43 


28.12 


0.48 


TXS 08284 


-193 


RL 


118 


2.577 


0.7x0.5 


44.51 


6.6 


0.56 


26.66 


0.61 



Notes. - All notations are identi cal to Table [T] and [2] R edshift (determined from the [O lll]A5007 A line), seeing and [O 
III] luminosity data are taken from lNesvadba et al.l l|2008t ). Radio fluxes are from the NVSS survey l|Condon et ahlfiggsh . 
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